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<^ [ Abstract 

We investigate hadron formation in deep inelastic lepton scattering on N, Kr and 
Xe nuclei in the kinematic regime of the HERMES experiment. The elementary 
electron-nucleon interaction is described within the event generator PYTHIA while 



a full coupled-channel treatment of the final state interactions is included by means 
of a BUU transport model. We find a good agreement with the measured charged 
hadron multiplicity ratio for TV and Kr targets by accounting for the deceler- 
ation and absorption of the primarily produced particles as well as for the creation 
of secondary hadrons in the final state interactions. 
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Hadron production in deep inelastic lepton-nucleus scattering (DIS) offers a 
promising tool to study the physics of hadronization (1). The reaction of the 
exchanged virtual photon (energy u, virtuality Q 2 ) with a bound nucleon leads 
to the production of several hadrons. While the primary production is deter- 
mined by the fragmentation function - in medium possibly different from that 
in vacuo - the number of ultimately observed hadrons and their energy distri- 
bution depends also on their rescattering in the surrounding nuclear medium. 
Consequently, the particle spectrum of a lepton-nucleus interaction will differ 
from that of a reaction on a free nucleon. In order to explore such attenu- 
ation effects the HERMES collaboration has investigated the energy v and 
fractional energy Zh = E^/v dependence of the charged hadron multiplicity 
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in DIS off iV (2) and Kr (0) nuclei. Here N^iz^u) represents the number of 
semi-inclusive hadrons in a given (z, v)-hm. and N e (v) the number of deep 
inelastically scattered leptons in the same z/-bin. It was suggested in Ref. (@), 
that a phenomenological description of the R% data can be achieved if the 
formation time, i.e. the time that elapses from the moment when the pho- 
ton strikes the nucleon until the reaction products have evolved to physical 
hadrons, is assumed to be proportional to (1 — z^)v in the target rest frame. 
This (1 — Zh)v dependence of the formation time r/ is compatible with the 
gluon-bremsstrahlung model of Ref. (Jl|). In the investigations of Ref. (j2j) any 
interaction of the reaction products with the remaining nucleus during this for- 
mation time has been neglected. After the formation time the hadrons could 
get absorbed according to their full hadronic cross section. Another interpre- 
tation of the observed R h M spectra - as being due to a combined effect of a 
rescaling of the quark fragmentation function in nuclei due to partial decon- 
finement as well as the absorption of the produced hadrons - has recently been 
given by the authors of Ref. (0). Furthermore, calculations based on a pQCD 
parton model (jD; 0) explain the attenuation observed in the multiplicity ratio 
solely by partonic multiple scattering and induced gluon radiation completely 
neglecting any hadronic final state interactions (FSI). It has already been 
pointed out by the authors of Ref. (J^y) that a shortcoming of the existing mod- 
els is the purely absorptive treatment of the FSI. We will avoid this problem 
by using a semi-classical coupled-channel transport model (0) which already 
has been employed to describe high energy photo- (8) and electropro duct ion 
off nuclei. 



We point out that the formation time also plays an important role in studies 
of ultra-relativistic heavy-ion reactions. For example, the observed quenching 
of high transverse momentum hadrons in Au + Au reactions relative to p + p 
collisions is often thought to be due to jet quenching in a quark gluon plasma. 
However, the attenuation of high vr hadrons might at least partly be due to 
hadronic rescattering processes (|lOt 111). 



Unfortunately, the average hadron formation time is not well known and the 
number of T/ ~ 1 fm/c as used commonly in the Bjorken estimate for the 
energy density |l2[ ) is nothing but an educated guess. The nonperturbative 
nature of this number - due to time scales of ~ 1 fm/c and hadronic size 
scales of 0.5-1 fm - excludes perturbative evaluation schemes; it is hard to 
calculate 7/ from first principles and formation times cannot be addressed in 
present lattice QCD simulations. One might expect that the rather successful 
string models (|13T ) shed some further light on this number, since the intrinsic 
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time scale tq for the q — q formation vertex can be related to the fragmentation 



function and string tension, respectively However, the actual parameters 
employed in current transport codes are not unique, with hadron formation 
times ranging from 0.3-2 fm/c (0E3), 

depending on the flavor, momentum 
and energy of the created hadrons. In fact, the rapidity and transverse mass 
spectra from relativistic nucleus-nucleus collisions are not very sensitive to the 



formation time t/ It is therefore essential to check whether these times 
are compatible with constraints extracted from reactions, where the collision 
geometry is much better under control. 

To this aim the attenuation of antiprotons produced in p + A reactions at 
the AGS energies of 12.3 GeV and 17.5 GeV has been investigated on various 
nuclear targets in Ref. (fl^i and a range of values for Tf = 0.4 — 0.8 fm/c has 
been extracted in comparison to the data from the E910 Collaboration (fill). 

In our present approach the lepton-nucleus interaction is split into two parts: 
1) In the first step the virtual photon is absorbed on a bound nucleon of the 
target; this interaction produces a bunch of particles that in step 2) are prop- 
agated within the transport model. Coherence length effects in the entrance 
channel, that give rise to nuclear shadowing, are taken into account as de- 
scribed in Ref. (0). The virtual photon-nucleon interaction itself is simulated 
by the Monte Carlo generator PYTHIA v6.2 (19) which well reproduces the 
experimental data of Ref. f)20r ) on a hydrogen target. Depending on whether 
the photon interacts directly with the nucleon or via one of its hadronic fluctu- 
ations (p°, u, J/ty or a perturbative qq fluctuation) the reaction leads to the 
excitation of one or more hadronic strings which in our approach are assumed 
to fragment very rapidly into colorless prehadrons. 

The time that the reaction products need to evolve to physical hadrons, i.e., 
the production time t p of the prehadrons plus the time needed to build up 
the hadronic wave function, we denote as formation time tf in line with the 
convention in transport models. For simplicity we assume that the formation 
time is a constant r/ in the rest frame of each hadron and that it does not 
depend on the particle species. Due to time dilatation the formation time tf 
in the laboratory frame is then proportional to the particle's energy 

+ ZhV <n\ 
t f = nf. T f= Tf. (2) 



The size of t/ can be estimated by the time that the constituents of the hadrons 
need to travel a distance of a typical hadronic radius (0.5-0.8 fm). 

To illustrate the situation shortly after the photon nucleon reaction Fig. 1 
shows the excitation and fragmentation of a hadronic string in a deep inelastic 
scattering process. For simplicity we do not show any gluon bremsstrahlung of 
the struck quark in this figure. Note, however, that the possibility of such final 
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state gluon radiation plus subsequent qq splitting is included in the PYTHIA 
part of our model and leads to the creation of additional strings (fljf). It has 
been emphasized in Ref. (jjlh that the string propagating through the nucleus 
is a rather short (white) object of length 1 fm since the slow end of the 
string is accelerated very fast. When the primary string fragments - due to the 
creation of qq pairs from the vacuum - new colorless prehadrons are produced, 
which we propagate in space-time. 

As discussed in Ref. (J22I ) the production time t p of these prehadrons has to 
be distinguished from the total formation time tf of the final hadrons which 
is dilated according to Eq. (2). It has been confirmed (|23h within the Lund 
model that the production time t p vanishes for Zh — > and z^ — > 1. In the 
present numerical realization of our model we first approximate this behavior 
by setting the production time t p to zero for all prehadrons, but will also 
discuss the effect of a finite production time at the end of this study. 



Right after the photon nucleon interaction the primary string should interact 
with a hadronic cross section because its transverse size is essentially that 
of the original nucleon. Motivated by the constituent quark model we as- 
sume that this hadronic cross section is shared by the quark/diquark at the 
string ends and after the fragmentation by the so called leading prehadrons 
that contain this quark or diquark. Our PYTHIA simulations show that in 
most cases the prehadrons with Zh ~ 1 are such leading hadrons since they 
contain constituents (valence- or sea-quarks) from the target nucleon or the 
hadronic component of the photon. They can therefore interact directly after 
the photon-nucleon interaction with a constant effective cross section which 
we denote as cri ea d- The cross sections of the other prehadrons, that solely 
contain quarks and antiquarks created from the vacuum, emerge at interme- 
diate Zh. They are assumed to be non-interacting until tf. This assures that 
the summed cross section of the complete final state right after the photon- 
nucleon interaction is approximately that of the original nucleon. Each time 
when a new hadron has formed, the summed cross section rises just like in 
the approach of Ref. (21). After the hadron formation time tf by definition all 
hadrons interact with their full hadronic cross section a^. Note, that our con- 
cept of leading hadrons is in accordance with those of other transport models 
for high energy reactions ([11 : 



Since the lighter (intermediate Zh) hadrons have large formation times in the 
target frame (see Eq. (2)) they may escape the nucleus without being attenu- 
ated if they are non- leading. However, many (« 2/3) of the observed hadrons 
with intermediate are not directly produced in the string fragmentation 
but stem from decays of the much heavier vector mesons with correspond- 
ingly shorter formation times. These vector mesons may therefore form inside 
the nuclear volume and thus be subjected to FSI. The effect of the FSI, fi- 
nally, will depend dominantly on the nuclear geometry, i.e. the size of the 
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target nucleus. 

In our present study the FSI are described by a coupled-channel transport 
model based on the Boltzmann-Uehling-Uhlenbeck (BUU) equation which de- 
scribes the time evolution of the phase space density fi{r,p,t) of particles of 
type % that can interact via binary reactions. These particles involve nucleons, 
baryonic resonances and mesons (n, rj, p, K, ...) that are produced either in the 
primary reaction or during the FSI. In this work we also have to account for 
the prehadrons emerging from the string fragmentation. For a particle species 
i the BUU equation takes the form 

(0 OH OH \ 
m + w " ~d?w) hm t] = ~ fu /m] ' (3) 



where the Hamilton function H includes a position and momentum dependent 
mean-field potential for baryons. The collision integral on the right hand side 
accounts for the creation and annihilation of particles of type i in a collision 
as well as elastic scattering. The transition rates are determined from the 
particular (vacuum) cross sections. For fermions Pauli blocking is taken into 
account in I co u via blocking factors. The prehadrons are treated like ordinary 
hadrons except for their modified interaction cross section and the fact that 
they are not allowed to decay during the formation time. The BUU equations 
of each particle species i are coupled via the collision integral and the mean 
field in case of baryons. The resulting system of coupled differential-integral 
equations is solved via a test particle ansatz for the phase space density. For 
further details of the transport model we refer the reader to Ref. (0). 

Most of the hadronic FSI happen with invariant energies > 2.2 GeV and 
are described within the Lund string formation and decay schem e (1131) a s also 
implemented in the transport approaches (0;S0) as well as ([lit 1 171: 124). The 
important difference between a purely absorptive treatment of the FSI and 
the coupled-channel description provided by the BUU model is that in an in- 
teraction with a nucleon a hadron might not only be absorbed or recreated 
but also be decelerated in an elastic or inelastic collision. Furthermore, it may 
in addition produce several low energy particles. In the case of electroproduc- 
tion of hadrons these interactions leads to a redistribution of strength from 
the high Zh part of the hadron energy spectrum to lower values of the energy 
fraction z^. 

In our calculation we employ all kinematic cuts of the HERMES experiment as 
well as the geometrical cuts of the detector. In actual numbers: we require for 
the Bjorken scaling variable x = 2 m N v > ^-06 (with denoting the nucleon 
mass), for the photon virtuality Q 2 > 1 GeV 2 and for the energy fraction 
of the virtual photon y = u/E heam < 0.85. In addition, the PYTHIA model 
introduces a lower cut in the invariant mass of the photon-nucleon system at 
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W = 4 GeV that is above the experimental constraint W > 2 GeV. This limits 
our calculations to minimal photon energies of v min = 8.6 GeV as compared 
to v m in — 7 GeV in the HERMES experiment and leads to a suppression of 
high Q 2 events at energies below v « 15 GeV. 

In Fig. 2 we present the average values of Q 2 and v (zh) for our simulated event 
samples on N and Kr as a function of Zh (y) in comparison to the experimental 
quantities (J26|). In order to compare with the v dependence of (Q 2 ) and (z h ) 
we only account for hadrons with Zh > 0.2 as in the HERMES experiment. 
For both the N and the Kr target the average kinematical variables are well 
reproduced by our model. The underestimation of (Q 2 ) at low photon energies 
v < 12 GeV is due to the PYTHIA cut in W > 4 GeV which suppresses higher 
values of Q 2 at low photon energies v. This is also the reason why the average 
value of Q 2 in the Zh spectrum comes out slightly too low within our model. 

We proceed with a discussion of the actual results of our PYTHIA + BUU 
simulations. In Fig. 3 we show the calculated multiplicity ratio for iV 
and Kr targets using an 'estimated' formation time of 0.5 fm/c and different 
values for the leading prehadron cross section o"i ea d. The data have been taken 
from Refs. (0; 0). Since the particles with Zh close to 1 are predominantly 
leading hadrons we can use the high Zh part in the fractional energy spectrum 
to obtain information on o"i ea d- The data for both nuclei indicate that awi 
has to be in the range 0-0.5 Oh with ah (h = vr* 1 , K ± ,p, . . .) taken from ([271). 
For the heavier nucleus, Kr, we clearly underestimate the hadron attenuation 
with diead =0 since most of the particles, especially those with large energies, 
escape the nucleus due to time dilatation. If one wants to describe the strong 
attenuation of hadrons at large fractional energy Zh without any prehadronic 
interactions one would need an unphysical short formation time r/ < 0.1 
fm/c. This, however, is ruled out by the measured v dependence of R\ { since a 
vanishing formation time leads to a multiplicity ratio R h M which is considerably 
too low (see dash-dotted and dash-dot-dotted curves in Fig. 4). 

The dotted line in Fig. 3 shows the result of a calculation where the leading 
prehadrons interact with the full hadronic cross section a^', this leads to a too 
strong attenuation of charged hadrons. A good agreement with the data is 
achieved for <7i ea d = 0.330"^ during the formation time Tf as can be seen in Fig. 
4. We note, that this value for o"i ca d represents an average value over time from 
the virtual photon-nucleon interaction to the actual hadron formation time. 
For a detailed investigation we refer the reader to a forthcoming study (j28h . 

In Fig. 4 we investigate the influence of different formation times Tf on R\j 
using the effective cross section o"i ea d = 0.33(7^. We find, that formation times 
Tf > 0.3 fm/c are needed to describe the experimental data with little sen- 
sitivity to higher values. This is compatible to the range of values extracted 



from the antiproton attenuation studies in Ref. (17). The steep rise of R^ 



6 



for Zh < 0.2 is caused by the energy loss and the production of low energy 
secondary particles in elastic and inelastic FSI. We mention that in models 
dealing with purely absorptive FSI a ratio much larger than one can only be 
explained by a drastic change of the fragmentation function. 

We note that the photon energy dependence of the ratio R\j for the Kr target 
is less well reproduced for energies below 14 GeV. The reason for this can be 
traced back to the W > 4 GeV cut of the PYTHIA model. Our simulations 
show that the number of leading hadrons decreases with Q 2 for fixed energy, 
but the W> 4 GeV cut discards larger values of Q 2 at energies v < 15 GeV. For 
higher values of Q 2 the importance of DIS events rises as compared to events, 
where the photon interacts via a vector meson fluctuation (VMD events). In 
the latter case one has initially 5 constituent quarks and antiquarks and thus 
gets more leading prehadrons than in a DIS event. This leads to a mismatch of 
leading and nonleading prehadrons in the region Zh < 0.4. Since the number 
of prehadrons created in an electron-nucleon collision decreases exponentially 
with increasing z^, this region contributes dominantly to the z^ integrated v 
spectrum. This deficiency of the present model will be cured in a more detailed 
upcoming work (|28h. 

In the pQCD parton model of Ref. (0) the multiple parton scattering leads 
to a modification of the fragmentation function and predicts a hadron atten- 
uation ~ A 2 / 3 . In our approach, however, the fragmentation is assumed to be 
decided on time scales of the nucleon dimension itself such that only the 'free' 
fragmentation function enters. All attenuation effects then are attributed to 
FSI of the leading and secondary (pre-)hadrons. In order to distinguish ex- 
perimentally between the different concepts, it is thus important to get the 
scaling with target mass A. To this aim Fig. 4 also shows predictions for a 
Xe target. In accordance with the authors of Ref. (0) we predict only a small 
change in the multiplicity spectra compared to the Kr target such that the 
scaling exponent is lower than 2/3. 

Finally, we discuss the effect of a finite production time t p (in the lab frame) 
for prehadrons which we adopt, for consistency, from the Lund model (j23h: 



Here k k 1 GeVfm -1 denotes the string tension. In the following we assume 
no interaction before t p and the full hadronic cross section for prehadrons for 
t > t p , i.e., there is no dependence on the formation time tf anymore. The 
dashed line in Fig. 5 represents the result of a Glauber-like treatment of the 
FSI where every time a prehadron interacts with another particle it is removed 
from the outgoing channel. As the authors of Ref. (0) we get a good description 
of the z-dependence of the multiplicity ratio. However, the v spectrum is not 
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attenuated strongly enough at the higher v (r.h.s. of Fig. 5). The solid curves 
show the effect of the coupled channels, i.e., a particle is not only absorbed in a 
collision but produces a bunch of low energy particles, thereby shifting strength 
to the low z part of the spectrum and thereby underestimating the attenuation 
at low z. Similarly, the attenuation is also too weak in the //-spectrum. Since an 
additional formation time with reduced cross sections would further enhance 
these discrepancies, we conclude that the data cannot be described with the 
production time (4). For a detailed investigation of finite production times 
and alternative models we refer the reader to a forthcoming work (fiih. 

In summary, we have shown that one can describe the experimental data of 
the HERMES collaboration for hadron attenuation on nuclei without invok- 
ing any changes in the fragmentation function due to gluon radiation. In our 
dynamical studies, that include the most relevant FSI, we employ only the 
'free' fragmentation function on a nucleon and attribute the hadron atten- 
uation to the deceleration of the produced (pre-)hadrons due to FSI in the 
surrounding medium. We find that in particular the z-dependence of R\ { is 
very sensitive to the interaction cross section of leading prehadrons and can 
be used to determine cri ea d- The interaction of the leading prehadrons during 
the formation time could be interpreted as an in-medium change of the frag- 
mentation function, which however could not be given in a closed form. The 
extracted average hadron formation times of tj > 0.3 fm/c are compatible 
with the analysis of antiproton attenuation in p + A reactions at AGS energies 
|l7h. In an upcoming work we will investigate in detail the spectra for differ- 
ent particle species (7r ± , K ,p,p) to examine, if the formation times of mesons 
and antibaryons are about equal. In addition we will improve our model to 
describe the primary photon-nucleon reaction below the PYTHIA threshold 
of W > 4 GeV. 

The authors acknowledge valuable discussions with N. Bianchi, E. Garutti, C. 
Greiner, V. Muccifora and G. Van Steenhoven. This work was supported by 
DFG. 
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Fig. 1. Illustration of an electron nucleus interaction: The virtual photon 7* excites 
a hadronic string by hitting a quark q inside a bound nucleon. In our example the 
string between the struck quark q and diquark qq fragments due to the creation of 
two quark-antiquark pairs. One of the antiquarks combines with the struck quark to 
form a 'leading' pre-meson m, one of the created quarks combines with the diquark 
to form a 'leading' pre-baryon b. The remaining partons combine to a pre-meson 
w! that, depending on the mass of the meson, might leave the nucleus before it 
hadronizes (see Eq. (2)). Note that in our approach the actual production time 
t p of the non-leading prehadrons has no effect on our results since we neglect any 
interaction until tf. See text for details. 
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Fig. 2. Model predictions for the average values of the kinematic variables in charged 
hadron production in comparison with the experimental numbers at HERMES. For 
the calculation we used the formation time tj = 0.5 fm/c and a leading prehadron 
cross section 0"i ea d = 0.33(7^. Left: (v) and ( Q 2 ) as a function of Zh compared to the 
experimental values for N and Kr targets (|2fil ). Right: Same for (zh) and (Q 2 ) as a 
function of v. 
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Fig. 3. Calculated multiplicity ratios (Eq. (1)) of charged hadrons for N and Kr 
targets for fixed formation time t/ = 0.5 fm/c and different values of the leading 
prehadron cross section: o"i ea d = 0., i.e. without any prehadronic interaction (solid 
line), criead = 0.5(7/! (dashed line) and ai ea d = ah (dotted line). The data for the 
Nitrogen target have been taken from Ref. (J2) while the Krypton data stem from 
Ref. 0). 
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Fig. 4. Calculated multiplicity ratios of charged hadrons for N, Kr and Xe targets 
for a fixed leading prehadron cross section <ri ea d = 0.33c/! and different values of the 
formation time from ry = to 1.5 fm/c. The data are the same as in Fig. 3. 
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Fig. 5. Calculated multiplicity ratios of charged hadrons for a Kr target assuming 
the finite production time t p given by Eq. (4). For t > t p the prehadrons interact 
with their full hadronic cross section. The dashed line shows the result of a simu- 
lation with a purely absorptive treatment of the FSI. The solid line represents the 
coupled-channel result. 
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